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ABSTRACT 

in ■ 

o : 

In binary radio pulsar system J0737— 3039, slow pulsar B shows orbital modulations 
^ ■ of intensity, being especially bright at two short orbital phases. We propose that these 

modulations are due to distortion of pulsar B magnetosphere by pulsar A wind which 
produces orbital phase-dependent changes of the direction along which radio waves are 
\ emitted. In our model, pulsar B is intrinsically bright at all times but its radiation 

beam misses the Earth at most orbital phases. We employ a simple model of distorted 
B magnetosphere using stretching transformations of Euler potentials of dipolar fields. 
^r) ' To fit observations we use parameters of pulsar B derived from modeling of A eclipses 

(Lyutikov & Thompson 2005). The model reproduces two bright regions approximately 
■ at the observed orbital phases, explains variations of the pulse shape between them and 

. regular timing residuals within each emission window. It also makes predictions for 

, timing properties and secular variations of pulsar B profiles. 

Or 

Pi 1. Introduction 

Double pulsar system PSR J0737— 3039A/B contains a recycled 22.7 ms pulsar (A) in a 2.4 hr 
. orbit around a 2.77 s pulsar (B) (Burgay et al. 2003; Lyne et al. 2004). Emission of pulsar B is 

strongly dependent on orbital phase. It is especially bright at two windows, each lasting for about 
30°. One of the emission windows appears near superior conjunction (when pulsar B is between 
pulsar A and observer), and another approximately 70° before it. Pulse profiles have different 
shapes in the two windows. In addition, Ransom et al. (2004a) have detected regular, orbital 
phase-dependent drift of emission arrival times by as much as 20 ms. 

Previously, Jenet &: Ransom (2004) suggested that emission of B is initiated by 7-ray emission 
from A. This model seems to be inconsistent with the evolution of the A profile (Manchesterei al. 
2005). Zhang & Loeb (2004) suggested that B emission is triggered by particles from pulsar A wind 
reaching deep into magnetosphere. This model is incorrect since the authors neglected magnetic 
bottling effect, which would reflect most pulsar A wind particles high above in B magnetosphere 
(see Lyutikov & Thompson 2005, for discussion of particle dynamics inside B magnetosphere) 

In this paper we explore an alternative possibility that orbital brighting of B is due to dis- 
tortions of B magnetosphere by pulsar A wind. We show that due to orbital phase-dependent 
distortions of B magnetosphere, the polar magnetic field lines, along which emission is presumably 
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generated, may be "pushed" in the direction of an observer at particular orbital phases, while at 
other moments radiating beam of B misses the observer. 

2. Distortion of pulsar B magnetosphere 

Magnetosphere of pulsar B is truncated, if compared to magnetosphere of an isolated pulsar of 
the same spin, by the relativistic wind flowing outward from pulsar A. The size of magnetosphere 
is R m ~ 4 x 10 9 cm (Lyutikov 2004; Arons et al. 2004; Lyutikov & Thompson 2005), which 
is 3 times smaller than the light cylinder radius Rlq = = 1-3 x 10 10 cm (U is angular 
frequency of pulsar B rotation, c is velocity of light). At intermediate distances, Rns < r < 
(Rlc, Rm), neutron star magnetospheres are well approximated by dipolar structure. This has 
been a longstanding assumption in pulsar theory (Goldreich & Julian 1969), and has recently been 
confirmed by modeling of pulsar A eclipses (Lyutikov & Thompson 2005). 

A small degree of distortion of magnetosphere from the dipolar shape is expected at interme- 
diate distances due to several types of electrical currents. First, distortions are due to confining 
Chapman-Ferraro currents (Chapman & Ferraro 1930) flowing in the magnetopause (a region of 
shocked pulsar A wind around B magnetosphere). At the emission radius, R em ~ 1 — 5 x 10 8 
cm (see below), fractional distortions due to Chapman-Ferraro currents are expected to be small 
~ (R em /R m ) 3 ~ .001 — .02. In addition to confinement of magnetosphere on the side facing pulsar A 
("dayside"), on the side opposite to pulsar A ( "nightside" ) magnetosphere of B extends to large dis- 
tances, somewhat similar to the Earth magnetosphere under the influence of Solar wind. Secondly, 
similar to isolated pulsars, there are conduction Goldreich-Julian currents, arising on the open field 
lines due to relativistic electromagnetic effects of rotation, and displacement currents, arising in 
oblique rotators due to temporal variations of electromagnetic fields. At the emission radius these 
currents produce a distortion of magnetic field of the order ~ (R em /RLc) 2 ~ 10~ 3 . Thirdly, there 
are internal currents flowing in the magnetosphere, like ring current, Birkland currents, and other 
types of currents. 

Clearly, the detailed structure of pulsar B magnetosphere is even more complicated than that 
of the Earth, but at distances from the star smaller than light cylinder and magnetospheric radii 
distortions from dipolar form are expected to be small. This should be the case in the emission 
generation region. In addition, since magnetospheric radius R m is several times smaller than light 
cylinder we can make a simplifying assumption that at each given moment the structure of the inner 
magnetosphere is determined by the instantaneous direction of the magnetic moment of B and the 
direction of line connecting two pulsars (which is, approximately, the direction of pulsar A wind at 
the position of B). Under this approximation we expect that at each moment the structure of the in- 
ner magnetosphere may be estimated using methods developed for non-relativistic, quasi-stationary 
magnetospheres of solar planets interacting with the solar wind. This interaction is complicated, 
depending on a number of both macroscopic (e.g. wind pressure, direction of magnetic field, dipole 
inclination) and especially microscopic (e.g. reconnection and diffusion rates) parameters. In what 
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follows we use experience with the solar wind - Earth magnetosphere interaction (e.g. Tsyganenko 
1990) as a guiding line in studying pulsar B magnetosphere. 

One of the principal issues here is how magnetopause currents respond to the dipole magnetic 
field of the central object. Two extreme possibilities are (i) complete screening of the dipole, so 
that no pulsar magnetic field lines penetrate into the wind and (ii) efficient reconnection so that 
most of the pulsar magnetic field lines that reach the magnetospheric boundary penetrate into 
the wind. In the case of the Earth magnetosphere, though reconnection plays an important role, 
on average inter-penetration of magnetic field is at most a 10% effect (eg. Stern 1987). (Rates 
of reconnection are strongly dependent on the direction of the solar wind magnetic field. On the 
dayside reconnection occurs most efficiently near the cusps, where polar magnetic field lines reach 
magnetopause.) Numerical modeling of interaction of relativistic pulsar A wind with pulsar B 
shows qualitatively similar results Arons et al. (2004). Thus, as a first approximation, we may 
assume that magnetopause currents screen out pulsar B magnetic field. 

An additional source of magnetic field distortion is the ring current generated by particles 
trapped inside magnetosphere. Modeling of pulsar A eclipses (Lyutikov & Thompson 2005) implies 
that high density, relativistic plasma (most likely composed of pairs) is present on closed field lines 
of pulsar B. In addition to bouncing between magnetic poles due to effect of magnetic bottling, 
trapped particles drift along magnetic equator. Viewed from the north magnetic pole, positrons 
drift clockwise, electrons - counterclockwise, producing a ring current, which modifies the structure 
of magnetosphere. We expect that effects of the ring current are negligible at the emission radius. 
A typical drift velocity of charge carries is Ud ~ c 2 7o/w£i? m ~ 3 x 10 4 cm/sec (here 70 ~ 10 
is a typical Lorentz factor of trapped particles, ujb is cyclotron frequency). For particle density 
n ~ X m nGJ,m (naj,m is Goldreich- Julian density at the magnetospheric radius R m = 4 x 10 9 cm and 
\ m is multiplicity factor at R m ), the current density is j ~ X m nGj, m Ude ~ -Ijcj, and total current is 
/ ~ jR 2 m . Resulting magnetic field is B ring ~ I/(R m c) ~ jR m /c ~ X m (QR m /c)(u d /c)B m ~ .03B m . 
Deep inside magnetosphere, the magnetic field of the ring current is nearly constant, but the dipole 
field increases, so that at R em ~ 10 8 cm B r i ng /Bd ~ 10~ 5 . (Qualitatively, the magnetic field of the 
ring current would become comparable to the dipole field at i? m when energy density of trapped 
plasma is of the order of magnetic field energy density.) 

There are other types of currents that can modify field structure like Birkland and tail cur- 
rents (e.g. Tsyganenko 1990). Their influence on the structure of magnetosphere at intermediate 
distances is expected to be small and we neglect them here. Thus, we assume that the only cur- 
rents contributing to the distortion of magnetosphere are magnetopause Chapman- Ferraro currents 
which screen pulsar B field. This simplification allows us to use models developed for the Earth 
magnetosphere in order to estimate distortions of pulsar B magnetic field. 
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3. Modeling distorted magnetosphere of B 

There is extensive literature on modeling of the Earth magnetic field (e.g. Tsyganenko 1990). 
A number of analytical methods have been developed. For our purposes, we do not need to calculate 
the full structure of magnetosphere, but only to estimate variations in the position of polar magnetic 
field lines, where emission of B is presumably generated. For this purpose we employ the method of 
distortion transformation of Euler potentials (Stern 1994; Voigt 1981). A major advantage of the 
stretching model of magnetosphere is that it reproduces fairly well the structure of a tilted dipole 
(Stern 1994). 

Magnetic field can be described by two Euler potentials a and (3 (sometimes called Clebsch 
potentials): 

B = Vq x V/3 (1) 

so that magnetic field line is defined by an intersection of surfaces with constant a and f). Magne- 
tosphere of B enshrouded by magnetopause resembles a dipole field compressed on the dayside and 
stretched out on the nightside. The structure of the nightside magnetosphere can be approximated 
by stretching transformations of the Euler potentials a and (3. 

Let us choose a system of Cartesian coordinates in the tail of B magnetosphere, centered on 
pulsar B so that z' axis is along the line connecting pulsar A and B and axis x' is in the fx — z' 
plane, where fl is magnetic moment of B, see Fig. 1. Let the magnetic dipole be inclined at angle 
9^ to axis z' . The undistorted dipole Euler potentials are 

(-x'smO^ + z'cose^ + y' 2 
ao = RnsV— 



yV 2 + y' 2 + 



z 

i 
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A) = Rns arctan 17 — — (2) 

—x' sin 6fj, + z' cos 9^ 

Stretching transformations along z 1 axis are defined by potential f(z'), so that new Euler potentials 
are expressed in terms of dipolar ones: {a, f3} = {olq(z' = f(z')),i%(z = f(z'))}. A degree of 
stretching depends on f'(z'). In modeling of the Earth magnetosphere function f(z') is chosen to 
fit satellite data. Since we are interested in the distortions at one particular location (assuming that 
radio emission is generated in a narrow range of radii), we choose f'(z') = C =const < 1, similar 
to Voigt (1981) model of Earth magnetosphere. Thus, C measures the distortion of magnetosphere 
at the emission radius. 

Substituting z' — ► Cz' in Eq. (2) we find stretched magnetic fields in coordinates {x',y',z'} 

B = % {SCx'z cos ^ + (2CV 2 - y' 2 - z' 2 ) sin M , 
r' 

3y'C(z' cos Of,, + Cx' sin 0^, 

C (3Cx'z f sin 0^ + (2z f2 - C 2 x' 2 - y' 2 ) cos 0„) } (3) 
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where r' = \JC 2 x' 2 + y' 2 + z' 2 . Integrating along a magnetic field line in the y' = plane we find 
equation for magnetic surfaces: 

r (C cos 6* M sin — sin ^ cos 0') 2 
r (cos 2 6>' + C 2 sin 2 6»') 3/2 

where r = yj x' 2 + z' 2 and 6*' = arcsinx'/r' are polar coordinates aligned with z' and r is a 
parameter related to maximum extension of a field line. From Eq. (4) we find that polar field lines 
in the tail of stretched magnetosphere are defined by 

tan P = — tan 6> M (5) 

This provides an estimate of the deviation of polar field lines from the direction of magnetic dipole. 
Qualitatively, polar field lines are pushed toward z' axis. The method of field line stretching is only 
approximate and has limited applications. Its main drawback is that it offsets force balance, so that 
there is a non-vanishing Lorentz force in the new configuration. In addition, since the stretching 
method has been devised for magnetotail, it's not clear how well it reproduces a structure of the 
inner magnetosphere (in original Voigt (1981) model the stretching method is applied to tailward 
distances larger than approximately half the stand-off distance) . 



4. Orbital modulation of B 

Let us introduce another Cartesian system of coordinates x, y, z centered on pulsar B, so that 
its orbital plane lies in the x — y plane (see Fig. 1). The spin axis of pulsar B is inclined at an 
angle 9^ to the orbital normal, and at angle (ftn with respect to the x — z plane. The magnetic 
moment of pulsar B has a magnitude \i, is inclined at an angle \ with respect to ft and executes 
a circular motion with phase 4> ro t = fit, so that 4> ro t = corresponds to the magnetic moment in 
the Q — x plane. At a given orbital position, the unit vector along the direction of pulsar A wind 
is approximately \ w = {cos 4>, sin <p, 0}. 

In the observer frame the components of the unit vector p,(t) along instantaneous magnetic 
moment jl(t) are 

p, x = p%; fiy = cos 9n fiy +sin9np%; fn z = cos 0q $ - sin6>n fay. (6) 
where fi are coordinates in a system aligned with CI, 

Ax = sin X cos((f> rot ); fly = sin x sin(^ rot ) ; fi^ = cos X- (7) 

Using Eq. (5) we can find the direction of magnetic polar field lines in the distorted magneto- 
sphere: 

A - (1 - C) cos 9 fXu 



VI + (3 - AC + C 2 ) cos 2 



(8) 
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where cos 9^ = ■ \ w \ is the absolute value of cosine of the angle between the direction of the wind 
and the direction of the magnetic moment. 

To estimate an influence of the orbit-dependent magnetic field distortion on observed radio 
emission, we assume that emission is generated near the polar field line in a region with half 
opening angle of ~ 2 degrees, in accordance with the width of pulsar B emission profile. Thus, 
if the magnetic polar field line deviates from the line of sight by more that 2 degrees, we expect 
emission of B to be weak. The trajectory that magnetic polar field line makes on the sky depends 
on many parameters. To limit available phase space we use the results of modeling of pulsar A 
eclipse (Lyutikov & Thompson 2005), which imply that 9q ~ 60° and x ~ 75°. In total, we have 
to fit at least 6 parameters: (j)n,9n,x, magnetospheric radius, impact parameter and distortion 
coefficient C. Using results of Lyutikov & Thompson (2005), the two principal parameters that 
remain to be determined are <pn and C. For a given set of 0n,4>n,X an d C, the direction of the 
polar field line executes a non-circular curve on the sky. An observer will see emission when for 
some values of pulsar B rotation phase (p ro t and orbital phase </> or 6 the polar field line points within 
two degrees of the line of sight. 

Searching through parameter space we were trying to reproduce two emission windows located 
in the tail of B magnetosphere. After a number of trials, our best fit parameters are <j>n = —67.5° 
and C = 0.7. To illustrate the fit, in Fig. 2 we plot a value of cos o b = s p ■ x as a function of orbital 
phase 4> or b and rotational phase <f) ro t- Only points located above the line cos 2° = .9994 produce 
a pulse of radio emission. (We restrict ourselves to — ir/2 < 4>orb < 7r /2, corresponding to the 
tail pointing towards the observer.) Clearly, in the tail of magnetosphere the polar field lines are 
pointing toward an observer at two orbital phases separated by approximately 70°, see Fig. 4. One 
of the phases, nearly coincident with the conjunction, is located at orbital phase 4> or b = 1.7°. ..24° 
and the second is located at (p or b = 62°... 83°. Considering the simplicity of the model and the fact 
that we had to make a multi-parameter fit, the agreement with observations is impressive. In order 
to illustrate the effect of field distortion, in Fig. 3 we plot a trajectory of magnetic polar lines on 
the sky, folded over ir radians to show both poles. Emission is seen only at points approaching the 
direction to the observer within two degrees. 

Though we were able to obtains a satisfactory fit, we cannot guarantee that there is no other 
island in the phase space of 6 parameters that also satisfies these criteria. In addition, there are 
intrinsic ambiguities in the model, related to prograde versus retrograde rotation of pulsar B and 
to which magnetic pole is seen by the observer. 

5. Discussions and Predictions 

We have constructed a simple model of orbital variations of pulsar B emission. Small distortions 
of the inner magnetosphere of B by pulsar A wind change direction of the polar field lines, pushing 
radiative beam of B towards the line of sight at two orbital phases. The main implication of the 
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model is that B is always intrinsically bright. The model reproduces fairly well absolute orbital 
location of bright emission windows, their width ~ 20° and separation ~ 70°. It also naturally 
explains different profiles at two emission windows, since at different orbital phases the line of sight 
crosses the emission region along different paths. 

In this paper we considered only the nightside magnetosphere. The stretching method does 
not produce a realistic structure of the dayside (Stern 1994). Qualitatively, we expect that on 
the dayside polar field lines will also be shifted from the direction of the magnetic pole and will 
be pushed out of the line of sight, producing a dip in the light curve of B close to the inferior 
conjunction. Since emission beam is very narrow and distortions can be considerable, it is fairly 
easy for the beam of B to miss an observer. A more detailed model of magnetosphere is deferred 
to a subsequent paper. 

Our suggestion that pulsar B is always intrinsically bright is consistent with its spin-down and 
radio energetics. First, assuming that an extension of the last open field line is determined by the 
size of magnetosphere and that typical current density flowing on the open field lines is of the order 
of the Goldreich-Julian current density, the total potential over open field lines (a quantity that is 
usually related to efficiency of pair production, (e.g. Arons k, Scharlemann 1979)) is independent 
of R m : <&tot ~ ^Lsd/^kc, where L$d = 1-6 x 10 30 erg/s is the spindown luminosity of B Lyne 
et al. (2004). This ranks it as the 20th smallest (but not exceptional) out of nearly 1500 pulsars 
with measured spindown luminosities (see www.atnf.csiro.au/research/pulsar/psrcat). Secondly, 
its peak luminosity of ~ 3 mJy at 820MHz (Ransom, priv. comm.) is typical for isolated pulsars 
with similar properties. 

There is a number of predictions of the model. First, at different orbital phases the center 
pulse corresponds to somewhat different rotational phases. Near <fi or i, = 1.7°, the center pulse is 
at 4> ro t ~ 11.5°, decreasing to 4> ro t ~ 8° at or {, = 24°, remaining the same at or b = 62.5° and 
increasing back to 4> rot ~ 11.5° at 4> or b = 83.1°, see Fig. 5. Thus, one expects a drift of the profile as 
a function of an orbital phase. Particular values for the drift angle and rate are strongly dependent 
on the precise parameters of the model, but the type of evolution is generic: during a visible phase 
the profile drifts approximately by its width. If averaged over a bright emission window, the drift 
of the emission phase may be interpreted as a large timing noise of B. 

A drift of the emission phase as a function of orbital position may have already been observed. 
Ransom et al. (2004a) reported a systematic change in arrival times of B pulses by 10-20 ms. 
This is consistent with the prediction of the model, since a change in arrive phase by 2 degrees of 
rotation phase of B corresponds to ~ 15 ms. 

In our model, the weak emission of B observed throughout the orbit has a different origin than 
the bright emission. For example, it can be generated in a much wider cone, akin to interpulse 
emission (bridges) observed in regular pulsars. 

Our second prediction is related to a secular evolution of pulsar B emission properties. The 
emission beam of B is fairly narrow, so that small changes in the orientation of rotation axis of B 
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may induce large apparent changes in the profile. One possibility for the change in the direction of 
the spin is a geodetic precession of B, which should happen on a relatively short time scale ~ 70 yrs 
(Lyne et al. 2004). The geodetic precession will affect mostly the angle 4>n. From our modeling 
we find that changes of (J)q by as little as ~ 1° strongly affect observed B profile (see Fig. 5). 
Thus, we expect that profile of B may change on a times scale of less than a year. According to 
the model, average profile width remains approximately constant. Changes of the orbital phases 
of emission are not accompanied by substantial changes in the emission phase. (Note that if the 
emission geometry is non-trivial, e.g. elliptic instead of circular, one does expect changes in the 
emission phase.) Since at different epochs line of sight passes through different emission regions 
one may expect variations in pulse intensity. Thus, using different slices taken at different epochs 
one can construct a detailed map of the emission region. This should prove a valuable method in 
constraining pulsar radio emission mechanisms. 

A longer evolution of the profile cannot be predicted unambiguously since we do not know 
the direction of the drift. Two possibilities include increasing or decreasing \4>q + 90° |, see Fig. 6. 
In one case, two emission regions get closer together merging in one, while in the other case they 
separate and a new one appears approximately at a mirror reflection of the first, at (j) or t, ~ —90°. 1 
We would like to stress again that exact details of the secular evolution are hard to predict using 
this very simple model. 

For the stretching coefficient C = 0.7, the relative deformation of the magnetosphere at the 
emission radius is ~ 30%. This is a fairly large distortion, which favors large emission altitudes, 
Rem > 10 8 cm. Large emission altitudes in isolated pulsars have been previously suggested by 
Lyutikov et al. (f998). A more precise modeling of B magnetic field may reduce required distortion 
and thus allow somewhat smaller R em - Still, near stellar surface distortions are expected to be tiny, 
so that the model requires high emission altitude. 

The success of this simple model is somewhat surprising, given the fact that in all we had to 
fit many parameters with a required precision of 2 degrees. Qualitatively, the reason for the success 
of this model, as well as that of Lyutikov & Thompson (2005), is that, to the first order, magnetic 
field is well approximated by dipolar structure. Given the simplicity of the model, some of the 
parameters may not be well determined: small variations in parameters may induce relatively large 
observed changes. 

A number of effects may complicate the picture. First of all, a non-trivial geometry of the 
emission region, e.g. elliptical instead of circular, may increase the quality of the fit. The fact that 
B profiles in the two emission windows are different implies that emission geometry is indeed non- 
circular. A double hump profile in one of the windows also points to a more complicated emission 
geometry. Secondly, non-spherical A wind will produce distortions dependent on orbital phase. 



1 According to the model, we are not likely to lose pulsar B in the coming year, yet the model is not sufficiently 
detailed to guarantee it. 
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Also, if reconnection between wind and magnetospheric field lines is important, the structure of 
magnetosphere may depend on the direction of the wind magnetic field . We plan to address these 
issues in a subsequent paper. Clearly, a detailed modeling of B magnetosphere is required to further 
finesse the model. 

When the paper has been mostly completed we learned the results of Burgay et al. (2005), 
who found secular changes in B profiles in general agreement with predictions of the model. 

I would like to thank Maura McLaughlin and Ingrid Stairs for numerous enlightening discus- 
sions and comments on the manuscript. I also thank Gerry Atkinson, Marta Burgay, Joeri van 
Leeuwen and Scott Ransom. 
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Fig. 1. — Geometry of the model. Observer is located at x — > oo. The spin axis of pulsar B is 
defined by angles 9q and 4>q. The direction of the magnetic moment \x makes angle x with the 
spin axis and rotates around 0, every 2.77 s. A line connecting two pulsars lies in the x-y plane 
and makes angle 4> or b with x axis. (Note that we define orbital phase with respect to the point of 
superior conjunction and not with respect to ascending node. Two definitions differ by 270°.) Tail 
coordinates x', y', z' are also shown with vector jl in the x' — z' plane. Angle 9^ is between z' and 
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Fig. 2. — Cosine of the angle 9 ob between the line of sight and the direction of the polar field 
line in a distorted magnetosphere as a function of rotational phase <f> ro t = fit and orbital phase 
— 7r/2 < or {, < it/2 measured from the superior conjunction. Parameters of the model are Oq = 60°, 
(/)q = 67.5°, x = 73.6°, C = 0.7. Emission is visible if O & < 2° corresponding to cos# ;, > .9994. 
This happens at two orbital phases 4> orb = 0.03.. .0.43 (1.7°. ..24°) and <j) orb = 1.09.. .1.45 (62°. .83°). 
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12 3 

9, angle with respect to orbital plane normal 

Fig. 3. — Trajectory of the polar field line on the sky for fa = ir/8 = 22.5°, C = .7. Black dashed 
line: undistorted dipole, green solid is at orbital phase (p or b = — ir/4 (misses the observer), blue line 
is at orbital phase cj) or b = and red line is at orbital phase 4> or b = 7r/4. The location of observer is 
denoted by the circle. Upper and lower sets of curves correspond to two magnetic poles. 
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Fig. 4. — Configuration of the system, after Lyne et al. (2004). Light shades segments indicate 
orbital phases where B emission is strongest, dark regions indicate the location of emission in the 
best fit model. 
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Fig. 5. — Dependence of the rotation phase of emission on the orbital position and evolution of 
the emission pattern due to changes of cf>Q. B radio emission is seen at Earth at two shaded areas. 
At different orbital phases peak of emission occurs at different rotation phases. Typical variations 
in emission phase are of the order of the profile width. The corresponding drift of arrival times 
for pulsar B is ~ 15 ms. Region in black is the current best fit (f>n = —67.5°, region on red is for 
<Pn = -68.5°. 
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Fig. 6. — Secular changes in the position of the emission regions due to geodetic drift: (a) = 
—68.5°, (b) (j>ci = —65° (current best fit value is <fin = —67.5°). In case (a) the orbital separation 
between the two bright phases increases, in case (b) it decreases, so that two phases merge in one. 
Except for </>q, parameters are the same as in Fig. 2. 



